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SUMMARY 

A comparison of the secondary flows computed by the viscous Kreskovsky- 
Brlley-McDonald code and the Invlscld Denton code with benchmark experimental 
data for a turning duct Is presented. The viscous code Is a fully parabolized 
space-marching Navler-Stokes solver while the 1 nvl sc Id code Is a time-marching 
Euler solver. The experimental data were collected by Taylor, Whltelaw, and 
Ylanneskls with a laser doppler veloclmeter system In a 90° turning duct of 
square cross-section. The agreement between the viscous and 1 nvl sc 1 d compu- 
tations was generally very good for the streamwlse primary velocity and the 
radial secondary velocity, except at the walls, where slip conditions were 
specified for the Invlscld code. The agreement between both the computations 
and the experimental data was not as close, especially at the 60,0° and 77.5° 
angular positions within the duct. This disagreement, was attributed to In- 
complete modelling of the vortex development near the suction surface. 


INTRODUCTION 

The continued emphasis on secondary flow problems In the design and 
analysis of turbomachinery requires three-dimensional flow prediction capable 
of accurately modelling secondary flows. Three-dimensional elliptic Navler- 
Stokes codes such as those developed by Hah (ref. 1) have been shown to predict 
cascade secondary flows fairly well. However, such codes require large 
amounts of computer storage and long computation times, thus precluding their 
use for routine design and analysis of turbomachinery geometries at the present 
time. This situation Is expected to change with further developments In com- 
puter architecture and algorithm optimization. Meanwhile, three-dimensional 
parabolized Navler-Stokes codes such as the Kreskovsky-Brl ley-McDonald code 
(ref. 2) and three-dimensional Euler codes such as the Denton code (ref. 3) 
can provide valuable Information about secondary flows while requiring less 
computer storage and shorter computation times than an elliptic Navler-Stokes 
solver. 

This paper compares the secondary flows predicted by the Invlscld Denton 
code and the viscous Kreskovsky-Brl ley-McDonald code with experimental data 
for a 90° turning duct of square cross-section taken by Taylor, Whltelaw, and 
Ylanneskls (ref. 4). The experiment was performed In order to provide bench- 
mark data for use In computer code assessment. While the turning duct Is a 
rather simple geometry compared to a turbomachine blade row, It allows much 
less cumbersome experimental technique to be used and fairly simple grid gen- 
eration. In addition, the secondary flows generated In a turning duct are 
less complex than those found In a blade row. The results presented In this 
paper consist of streamwlse primary and radial secondary velocity profiles and 
secondary flow vector plots at various angular positions In the duct. 


SYMBOLS 


1 /2 

De Dean number, De = (d/2r ) Re 

d hydraulic diameter 

h duct spanwlse height 

Re Reynolds number, Re = dV/v 

r radial coordlate 

r* normalized radial coordinate, r* * (r - r Q )/(r^ - r 0 ) 

r„ mean radius of curvature, r » (r, + r )/2 
c c ' 1 o 

r^ Inner radius of curvature (suction surface) 

r Q outer radius of curvature (pressure surface) 

V bulk mean velocity 

streamwlse velocity component 

U 

V r radial velocity component 

V z spanwlse velocity component 

w duct radial width 

z spanwlse coordinate 

z* mormallzed spanwlse coordinate, z* = z/h 

e streamwlse coordinate 

v kinematic viscosity 


METHOD 

Invlscld Computation 

The Invlscld code, described In reference 3, Is an explicit time-marching 
solution of the unsteady Euler equations In finite-volume form. A two-level 
multigrid technique Is used to reduce computation time while variable time 
steps and linear smoothing are employed to control Instability. The uniformly 
spaced grid used for the turning duct computations consisted of 21 points In 
the radial direction between the pressure and suction surfaces, 21 points In 
the spanwlse direction between the endwalls, and 75 points In the streamwlse 
direction. The grid Is shown In figure 1. 

The Invlscld code originally had provisions for Inlet total pressure shear 
In only the spanwlse direction; modifications were made to allow a full speci- 
fication of the Inlet total pressure field with shear In both transverse 
directions. This field was calculated from data presented In reference 4. 

Slip conditions at the walls were specified as the total pressure corresponding 
to 30 percent of the maximum freestream velocity for the Inlet plane. Pre- 
liminary computations with lower velocity slip conditions were subject to In- 
stability and would require large amounts of smoothing to obtain a converged 
solution with consequent numerical viscosity effects. The smoothing used In 
this computation was limited to the minimum required to prevent decoupling of 
the solution on adjacent grid points. 
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The uniform transverse grid spacing was considered to be adequate, since 
4 grid points were located within the boundary layer and more exacting resolu- 
tion at the walls would be meaningless, given the slip approximation required 
for an Invlscld computation. Since the 1 nvl sc 1 d code could not be run at In- 
compressible velocities for this case, the average velocity over the entire 
flow field was set at a Mach number of 0.3 to allow some density variation In 
the conservation variables while minimizing gross compressibility effects 
because the viscous computation was Incompressible and the experimental data 
was for a water tunnel. Typical CPU execution times for these computations 
using the 21 x 21 x 75 grid were 60 to 75 min on an IBM 3033 and 10 to 15 min 
on a CRAY-1 . 


Viscous Computation 

The results from the viscous computation were taken from reference 2. The 
viscous code Is based upon a primary-secondary velocity decomposition method 
which produces a set of approximate governing equations to correct an Invlscld 
potential flow solution for viscous effects, secondary flows, and pressure 
field distortion. Finite-difference formulations of the correction equations 
are solved as an Initial-value problem In space by an Implicit single-pass 
forward-marching technique. The method assumes that streamwlse diffusion Is 
negligible and that the bulk pressure correction Included In the streamwlse 
momentum equations Is uniform over each cross section; thus, It may be cate- 
gorized as a "fully parabolized Navler-Stokes" solver (ref. 5). The grid used 
consisted of 30 radial points, 20 spanwlse points between the endwall and the 
plane of symmetry and 25 streamwlse points between e = 0° and 0 = 77.5°. 

The transverse grid was locally refined near the walls, typical CPU execution 
times for this computation were approximately 4 min on a COC 7600, which would 
correspond roughly to 7 to 8 min on an IBM 3033 or 1 to 2 min on a CRAY-1. 
These times do not Include the execution time for the a priori potential 
solution. 


Experimental Data 

The experimental data were taken from reference 4, where laser doppler 
veloclmetry was used to measure primary and secondary flows In water through a 
90° turning duct of 40 mm x 40 mm square cross-section with a curvature ratio, 
r c /w, of 2.3. As a consequence of the uniform cross-sectional area, both 
curvature and pressure gradient effects were present In the duct. The Inlet 
boundary layers on the pressure and suction surfaces and on the endwalls were 
found to be 15 percent of the hydraulic diameter at a Reynolds number of 40 000 
and a Dean number of 18 700. Measurements were obtained at the Inlet and exit 
survey planes shown In figure 2 and at angular positions of 30. 0°, 60.0°, and 
77.5° In spanwlse Increments of 5 percent and radial Increments of 20 percent 
beginning at 10 percent away from the pressure surface. 


RESULTS AND DISCUSSION 
Qualitative Flow Features 

The secondary flow field computed by the Invlscld code Is shown In 
figures 3 to 5 for the 30.0°, 60.0°, and 77.5° positions respectively. 
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These secondary flow vectors are computed as vector deviations from the grid 
lines, but should approximate the true secondary flows as long as the grid 
lines approximate the true streamlines. The development of two symmetrical 
vortices about the midspan position Is clearly dlscernable at the 30.0° posi- 
tion. The secondary flow moves along the sidewalls of the duct from the pres- 
sure side to the suction side In the positive radial coordinate direction. 

The magnitude of the secondary flow Increases significantly as the flow moves 
downstream to the 60.0° and 77.5° positions. At the 60.0° position, shown In 
figure 4, the vortex cores have been Intensified and at the 77.5° position, 
shown In figure 5, the magnitude of the velocities and the extent of the 
vortex cores have been strongly developed. The secondary flow field at the 
duct outlet 100 mm downstream of the bend Is shown In figure 6, where the 
magnitude of the secondary velocities has diminished and the two vortices have 
become slightly asymmetrical. The asymmetry may be caused by grid skew that 
Is Inherent In the grid generation algorithm used with the Invlscld code. 

In order to make quantitative comparisons of the Invlscld and viscous 
computational results with the experimental data, plots of both streamwlse 
primary and radial secondary velocities were made at fixed distances from the 
pressure surface to the suction surface (r* = 0.1 to r* = 0.9) as a function 
of distance from the endwall to the midspan position (z* » 0.0 to z* = 0.5). 
These are discussed In the following sections. 


Primary Velocity Profiles 

The streamwlse primary velocity profiles from the Invlscld and viscous 
computations and the experimental data are shown In figures 7 to 9 for angular 
positions of 30.0°, 60.0°, and 77.5° within the duct. Plots are shown In each 
figure for constant values of r*, beginning near the pressure surface at 
r* = 0.1 and ending near the suction surface at r* = 0.9. The velocities 
are normalized by the bulk mean velocity, V. In general, the agreement between 
the Invlscld and viscous computations Is very good, except near the walls where 
slip conditions were required for the Invlscld code. The computations match 
the experimental data closely near the pressure surface for r* = 0.1 to 
r* = 0.5; however, near the suction surface for r* = 0.7 to r* = 0.9, the 
agreement Is not as close, especially at the 60.0° and 77.5° positions. This 
can be attributed to large distortion of the primary flow by strong vortex 
flow that Is modelled Incompletely by both the Invlscld and viscous computa- 
tions. The same trends are Indicated In figure 10, where the experimental 
data and the Invlscld computation for the duct outlet position are presented. 
Results from the viscous computation were not available since It was terminated 
at the 77.5° position In reference 2. 


Secondary Velocity Profiles 

The radial secondary velocity profiles from the Invlscld and viscous com- 
putations and the experimental data are shown In figures 11 to 13 for angular 
positions of 30.0°, 60.0°, and 77.5° wUhln the duct. The velocities are 
normalized by the bulk mean velocity, V. Both computations match the ex- 
perimental data near the pressure surface, but the agreement becomes progres- 
sively worse near the suction surface. The Invlscld computation overpredicts 
the radial velocities at the endwalls (z* = 0.0) due to the slip boundary con- 
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dltlons, This effect Is very pronounced from r* * 0.5 to r* * 0.9 for the 
60.0° position (fig, 12) and from r* a 0.7 to r* * 0.9 for the 77.5° 
position (fig. 13). The strong reversals In the experimental data at r* * 0.9 
and z* a 0.35 for the 60.0° and 77.5° positions clearly Indicate vortex 
formation that both the Invlscld and the viscous computations do not show. 

The disagreement between the Invlscld computation and the experimental 
data near the suction surface extends over most of the duct height as the vor- 
tex flow becomes more strongly developed at the 77.5° position. This behavior 
Indicates that the Invlscld computation failed to reproduce the essential 
radial secondary flow features as the vortlclty Increases. In the presence of 
milder vortlclty at the 30.0° position, much better agreement was achieved. 

At the duct outlet position, shown In figure 14, the agreement between 
the Invlscld computation and the experimental data Is good at r* = 0.1. At 
r* « 0.3 to r* = 0.9, the agreement deteriorates much more than at the 77.5° 
position. The large discrepancies can be attributed to viscous effects which 
were not Included Iri the Invlscld code and which became predominant over cur- 
vature effects In the straight portion of the duct. Results from the viscous 
computation were not available at the duct outlet position since It was ter- 
minated at the 77.5° position In reference 2. 


Overall Comparison 

The comparison of the Invlscld computation with the viscous computation 
for a turning duct geometry was undertaken to assess the capability of the 
Invlscld code relative to the viscous code to predict the secondary flows In a 
geometry for which benchmark data existed. The Invlscld computation required 
approximately an order of magnitude more CPU e>ecut1on time than the viscous 
computation (exclusive of the potential solution) for a solution that, because 
of the Inherent Invlscld approximation, did not predict the secondary flows 
quite as well as the viscous solution. This Is to be expected, since the In- 
vlscld code Is an Iterative solver of the Euler equations while the viscous 
code Is a single pass solver of correction equations to an a priori potential 
solution. Other factors that should be considered In the overall comparison 
of the two codes Include the amount of Input data preprocessing required and 
the generality of the Input geometry requirements. The viscous code requires 
a three-dimensional potential solution to obtain streamline coordinates and an 
Initial pressure field and the version In reference 2 Is currently limited to 
ducts of rectangular cross section. The Invlscld code utilizes an Input gen- 
erator (ref. 6) that Is an order of magnitude faster than the potential solver 
used for the viscous code. The Input generator produces the geometrical coor- 
dinates from a minimum of geometry and flow variable Information by using cubic 
spline fits for Interpolation and extrapolation. The Invlscld code can be 
used for stationary or rotating axial-flow, radial-flow, or mixed-flow turbo- 
machine geometries, although small cusps may be required at the blade leading 
and trailing edges to minimize grid slope discontinuities. 


SUMMARY OF RESULTS 

The primary and secondary flows computed by the viscous Kreskovsky- 
Brlley-McDonald code and the Invlscld Denton code were compared with benchmark 
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experimental data for a turning duct. The viscous code Is a fully parabolized 
space-marching Navler-Stokes solver while the invlscld code Is a time-marching 
Euler solver. The experimental data were collected by Taylor, Whltelaw, and 
Ylanneskls with a laser doppler veloclmeter system In a 90° turning duct of 
square cross-section. The agreement between the viscous and Invlscld computa- 
tions was generally very good for the streamwlse primary velocity, except at 
the walls, where slip condition?; were specified for the Invlscld code. Both 
computations failed to predict the distorted primary velocity profile shown by 
the experimental data near the suction surface at the 60,0° and 77.5° posi- 
tions. The Invlscld computation overpredicted the radial secondary velocity 
at the endwalls due to the imposed slip condition. Both computations failed 
to predict the strong reversals In the radial secondary velocity profile shown 
by the experimental data near the suction surface at the 60,0° and 77.5° posi- 
tions. The Invlscld compuatlon failed to match the experimental data at the 
duct outlet position due to the predominance of viscous effects over curvature 
effects In the straight portion of the duct. The comparison of the computa- 
tions with the benchmark data has shown that both the viscous code and the 
Invlscld code can provide useful Information about secondary flow phenomena If 
the Inherent limitations of each code relative to an elliptic Navler-Stokes 
code are recognized. 
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Figure 2. - Duct geometry. 
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Figure 8. - Streamwise primary velocity profiles at 60. 0° position. 
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Figure 11. - Radial secondary velocity profiles at 30. 0° position. 
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